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A sufficient condition for entanglement in two-mode continuous systems is constructed with the 
interference visibility and the uncertainty of the total particle number. The observables to be 
measured (particle numbers and particle number variances) are relatively easily accessible exper- 
imentally. The method may be used to detect entanglement in light fields or in Bose-Einstein 
condensates. In contrast to the standard approach based on entanglement witnesses, the condition 
is expressed in terms of an inequality which is nonlinear in expectation values. The condition is 
constructed using uncertainty relations with the particle number and the destruction operators. 
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I. INTRODUCTION 

In spite of considerable efforts, the separability of gen- 
eral mixed quantum states is still an open problem, even 
if the whole density matrix is known. The positivity of 
the partial transpose ,1| of the density matrix is a neces- 
sary condition for separability; however, it is a sufficient 
condition only for the 2x2 (two qubits) and 2x3 dimen- 
sional cases and for two-mode Gaussian states in contin- 
uous variable systems. For higher dimensions there exist 
entangled states Q with positive partial transpose and 
the separability problem is not fully solved. For bipar- 
tite low-rank density matrices an operational criterion 
for separability is presented in Ref. Q. A necessary 
and sufficient condition for entanglement of all bipartite 
Gaussian continuous states is also known Q. Entan- 
glement can also be detected by several other methods, 
e.g., through witness operators whose expectation value 
is negative only for (some) entangled states Q- 

In an experiment the density matrix is usually not 
known, only partial information is available on the quan- 
tum state. One can typically measure a few observables 
and still would like to detect some of the entangled states 
. Thus to find a criterion for entanglement with easily 
measurable observables is crucial for entanglement de- 
tection. There are only few such criteria in the literature 
U S Eil El for detecting entanglement in complicated 
situations such as in many-particle or continuous variable 
systems. The method described in Ref. [jj detects entan- 
glement among cold atoms having two internal degrees 
of freedom based on inequalities with the total angular 
momentum components. Ref. [SJ derives a criterion for 
the entanglement between two pairs of bosonic modes in 
terms of the total angular momentum and the particle 
number. This criterion detects entangled states that are 
close to singlet states of two large spins. 

A criterion for detecting entanglement between two 
modes is given in 0, ^ . Ref. @ presents a scenario 
where one just has to measure the second moments of x 
and p for both systems. For example, if the inequality 



is fulfilled, then the state is entangled. Here (AX) 2 — 
(X 2 ) — (X) 2 , and Xk and pt are canonical operators satis- 
fying [xk,Pk] — i- Note that the left hand side of Ineq. 
is quadratic in expectation values. However, by setting 
(xfc) and (pk) to zero by single-party unitary operations, 
only terms linear in expectation values remain and the 
criterion is equivalent to an entanglement witness. 

A generalization of criterion IJ) detects all entangled 
two- mode Gaussian states |9|,ll0|. However, in many ex- 
perimental situations non-Gaussian states are prepared. 
For example, if one has N photons and sends them 
through a beam splitter or if one has N atoms in some 
internal state and applies an appropriate laser pulse, the 
state will be 



l*) = 



1 



V2^V! 



(at+6t)^|o,0) = 




ri, N—n) 

Here a and b are annihilation operators which are defined 
according to xa = (a + a')/^/2. This entangled state 
is not detected by the previous criterion as it will be 
demonstrated later. 

In this paper, we will present a criterion which: (i) 
requires measuring quantities which are easily accessible 
experimentally and; (ii) detects entangled states in the 
vicinity of state The paper is organized as follows. In 
Sec. II. the entanglement criterion is derived. In Sec III. 
it is discussed what states are detected by this criterion 
and issues concerning its possible experimental applica- 
tions are also considered. A summary is given in Sec. 
IV. 



II. ENTANGLEMENT CRITERION 

In this section we will show that for all separable states, 
i.e. states that can be written as 



P = ^2PkPk ® Pk \ 

k 



(3) 



(A(x A + x B )) 2 + {A( PA - PB )) 2 < 2 



(1) 



the following expression involving the variances of the 
total particle number N := a^a + b^b and of the operator 
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(a — b) is bounded from below as 

j(A p A0 2 + l}{(A P (a - b)f + l} > ^ + 1 (4) 

where (A p A) 2 := (Atyl) p - | (A) \ 2 (note that A need 
not be Hermitian). 

The physical motivation for this criterion comes from 
the observation, made in the context of Bose-Einstein 
condensates in Ref. ^1 (also explained in Sec. Ill), that 
it is not possible to have a fixed total particle number — 
corresponding to (A p 7V) 2 = — and perfect interference 
— corresponding to (A p (a— b)) 2 =0 — at the same time, 
unless the system under consideration is in a highly non- 
classical, i.e. entangled, state. 

In Sec. II. A a simple separability criterion will be 
proved. In Sec. II. B this criterion will be generalized. 
Technical details are in Appendix A. Sec. II. C proves 
our main result, Ineq. 

A. Simple criterion 

In this subsection a simple separability criterion will be 
derived based on uncertainty relations for the two sub- 
systems. In order to understand the connection between 
the uncertainty relations and the necessary condition for 
separability, it is instructive first to review how the condi- 
tion (JTJ was obtained starting from the single-subsystem 
uncertainty relations 0, 

(A p x A/B ) 2 + {A pPA/B ) 2 > 1. (5) 

For a separable state of the form the sum of uncer- 
tainties of the EPR type operators Xa + xb and pa — Pb 
can be written as 

(A p (x A + x B )) 2 + (A p (p A - pb)) 2 
= I>fc{(A p ^) 2 + (A p ^ A ) 2 

k 

+ (\bxb) 2 + (A pfPB ) 2 | 

+ J2 Pk \ ^ x a+ x b)) 2 p a^ p b- {{x A +x B )) 2 p 

k 

+ ({pa -pb)) 2 p a® p b - ({pa -Pb)) 2 p I > 2. (6) 

In the equality it was exploited, that for a product state 
the uncertainty of an EPR type operator splits into the 
sum of the corresponding single system uncertainties as 
(A p a^ p b(x a + x B )) 2 = (A p ax a ) 2 + (A p bx b ) 2 . Based 
on the uncertainty relations @ for the individual sub- 
systems, A and B, the first sum in Ineq. I© is bounded 
from below by 2. Since the second sum is non-negative, 
the left hand side of Ineq. ijfJJl is also bounded by 2 which 
finishes the proof of criterion (JTJ- Thus for separable 



states the sum of the variances of xa + x b and pa — Pb 
has the same lower bound as the sum of the correspond- 
ing single system uncertainties. Note that this simple 
relationship holds because the right hand side of the the 
uncertainty relation JSJ) is a constant. On the other hand, 
this bound is not valid for non-separable states. In this 
case the sum of the two uncertainties can even be zero 
[9j since xa + x B and pa — Pb commute. 

After reviewing the previous example with the EPR 
operators, we will prove that all separable states fulfill 

S p := (A P A0 2 + (A p (a b)) 2 > f({N) p ), (7) 

where 

f{N) = sjN^\+^--2. (8) 

This will be necessary later to obtain our main result, 
Ineq. @. Here (A p iV) 2 is the variance of the total 
particle number in the two-mode system (i.e., in a two- 
mode electromagnetic field or Bose-Einstein condensate 
in a double- well potential) while (A p (a — b)) 2 is related 
to the variance of the phase difference between the two 
modes. 

The proof of J7| is the following. For a separable state 
of the form © , the sums of the two uncertainties in 
can be written as S p = S P; o + £ p ,i where 

S Pl0 = ^p fe ((A p .iV yl ) 2 + (A p .a) 2 
k ^ 

+ (A p BN B ) 2 + (A p Bb) 2 + N 2 -(N) 2 p }, (9) 

S P ,i = 5>|((a-6)) p J 2 -|((a-6)) p | 2 . (10) 

k 

Here N A ■= a f a, N B := tfb and N k := {a? a + tfb) Pk . 
Using the Cauchy-Schwarz inequality, one can show that 
S p ,i > 0, and thus S p > £ p ,o- 

Next, we will need the following uncertainty relation 
proved in Appendix A. 1-2 

R p := (A P N A ) 2 + (A p a) 2 > L((N A ) p ), (11) 

where 



L(N) = 




Obviously, the same inequality is true for subsystem B. 
Ineq. (|ll|l is an alternative of the number-phase uncer- 
tainty without the problem of defining an appropriate 
phase operator and the difficulties due to the 2tt periodic 
nature of the phase [l3| . 

In our case the bound in the uncertainty relation 111|) 
is not a constant, but a function of an operator expecta- 
tion value. Thus the method presented for the EPR type 
operators cannot be used, and careful analysis of the dif- 
ferent properties of the function L(N) must be done [l4| . 
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Using inequality and the fact that the bound i|12|) 
fulfills L(N-l) + L(N 2 ) > L{N l + N 2 ) + L(0) we obtain 

S P > ^>{l(/V,) + L(0) + N 2 - (N) 2 p j. (13) 

Now using the fact that L(Nk) + iV| is a concave function 
of TVfc, we obtain £ p > L(N) + L(0) which proves Ineq. 
0. 

Condition J7J) corresponds to a line on the (A p iV) 2 - 
(A p (a — b)) 2 plane (solid line in the inset of Fig. 0. All 
separable states belong to points above this line and fulfill 
Ineq. 0) . Points below this line correspond to entangled 
states only. 

B. Generalization 

We would like to find more entangled states in the 
(A P N) 2 - (A p (a — b)) 2 plane. In order to do that we 
generalize J7J as 

£ Pl10 := w(A p N) 2 + (1 - w){A p {a - b)f > f w ({N) p ) 

(14) 

where < w < 1 determines the relative weights of the 
two terms and f w (N) is defined at the end of this sub- 
section. Ineq. (|14|) correspond a region above a line with 
slope TZuJ (dashed lines in the inset of Fig. 0. These 
lines are the tangentials of the curve enclosing all sepa- 
rable states. The points below this curve correspond all 
to entangled states. 

In order to obtain the lower bound fw{(N) p ), we have 
to follow a procedure similar to what was presented in the 
previous subsection. For a separable state one obtains 

V p , w > Y,Pk\w(A p AN A ) 2 + (l-w)(A pA a) 2 

k ^ 

+ w(A p BN B ) 2 + (l-w)(A p sb) 2 

+ wN*-w(N) 2 \, (15) 

In Appendix A. 3 we prove the following uncertainty 
relation 

R p , w := w(A p N A ) 2 + {l-w){A p a) 2 > L w ((N A ) p ), (16) 
where 

L W (N) = S V^-^+^ + f -5 XN>N Ll 
I (N - N L )w(l - w) ii N < N L . 

(17) 

Here Nl = (1 — w)/Aw. Ineq. 1)16(1 is the generalization 
of Ulljl for unequal weights for the two variances. For 
N < Nl the function L W (N) is linear and the slope is 
determined in such a way that there is not an abrupt 
change in the derivative of L W {N) at N = N^. 

In order to get a lower bound for S p w using the un- 
certainty relation (|16fl . one has to follow similar steps as 




FIG. 1: Numerical test of the inequality (|1J for the two-mode 
separability problem for N=200. (solid) Boundary of the re- 
gion defined by Ineq. JIJ. All states below this line are en- 
tangled, (diamonds) Points corresponding to separable states 
found numerically. Inset: (solid) Boundary of the region de- 
fined by Ineq. 0; (dashed) Boundary of the region defined 
by Ineq. (O for w = 0.3 and 0.7. 



in Sec. II. A. Using the facts that L w {Nx) + L W (N 2 ) > 
L W {N X + N 2 )+ L w (0) is fulfilled and L W (N) + wN 2 is 
a concave function of N, the lower bound for S PjU) is 
obtained as f w (N) = L W (N) + L w (0). 

C. Proof of main result 

In this subsection Ineq. 0} will be obtained by deter- 
mining the curve which has the lines corresponding to 
different w's as its tangentials. The tangentials of a hy- 
perbola (y+co) = C/(x+cq) are given by wx+(l~w)y = 
2^/w(l — w)C — Co. One can reformulate (|14|> by replac- 
ing the right hand side by a slightly weaker lower bound 
which fits this form 



f w {N) = ^w{\-w){N+-)-l. (18) 

Hence the equation for a hyperbola on the (A p iV) 2 - 
(A p (a — b)) 2 plane can be obtained (solid curve in Fig. 
01. Ineq. (0J corresponds to points above this hyperbola. 
Any state which violates this inequality is necessarily en- 
tangled. 

III. DISCUSSION 

First, the tightness of the necessary condition for sep- 
arability Q should be verified. Numerical checks show 
that it is quite strong (see Fig. 0). The diamonds indicate 
product states of the form |0) ®]\E r ) found numerically. 
The state in the origin of Fig. giving zero for both 
variances in Ineq. 10), is state 10) as can be shown as 
follows. Eigenstates of TV with TV particles have the form 
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l*> = J2 c n\n,N - n). The state (a - has N - 1 

particles, thus (a — b)\^) = A|^) is possible only if the 
eigenvalue A = 0. A state for which (a — b)\^f) = 0, has 
to fulfill c n+ i\/n + 1 = c n \/N — n. This determines the 
state uniquely as (0). 

Our method detects entangled states in the proxim- 
ity of state J2J on the (AN) 2 - (A (a - b)) 2 plane as 
shown in Fig. ^ (In this section we will omit the p in- 
dex.) Other interesting states on the (AN) 2 - (A(a—b)) 2 
plane: A separable state having (AN) 2 = is the con- 
vex combination of products of Fock states \rik) \N — rife). 
For these (A(a — b)) 2 = N. Separable states with per- 
fect destructive interference between the modes having 
(A(a — b)) 2 = are the convex combination of products 
of coherent states of the form afe)|afc + c) where c is 
a constant common for all product subensembles. For 
these states (AN) 2 > N 

According to our criterion, to detect entanglement in 
an experiment, the variances of N and (a — b) should 
be measured. A simpler scenario is to measure the vari- 
ance of N and ((at - &t)( a _ b)). Since (A(a - b)f < 
((at — b')(a — b)) the J7J inequality can be used for this 
case after replacing (A(a — bj) 2 by ((at — b' f )(a — &)). 
The latter is just twice the particle number in the mode 
b'=(a-b)/V2. 

The generation of state (J2J) is never perfect, thus the 
system is in a mixed state. Our method makes it possible 
to detect entanglement even in this case. If (b'^b 1 ) w 
(AN) 2 the maximum particle number variance for a state 
to be detected is (AN) 2 cx VN which is much smaller 
than for coherent states. On the other hand, for perfect 
destructive interference when {b''b') ~ the maximal 
variance is (AN) 2 cx N. 

Equ. describes the quantum state of a Bose- 

Einstein condensate of atoms, if the a and b modes cor- 
respond to the two halves of the condensate I n 
this case (at + 6+) creates a particle in state |\&) and @ 
describes a product of single particle states of the form 
|*) ® |*) ® • • • ® |*). Although it is a product state 
from the point of view of the individual particles, in the 
n,m) basis it is clearly entangled. In order to detect 
entanglement, one needs to measure the variance of the 
total particle number and the particle number in one of 
the new modes after the two halves of the condensates 
interfere fl5lll6l |. 

The condensate can be "split" into two modes, realiz- 
ing state and then reunited for detection in a Mach- 
Zehnder type interferometer [l6|. The state J2J) can also 
be obtained in a Bose-Einstein condensate of two level 
atoms, by preparing the atoms in the same internal state 
and then applying a n/2 laser pulse. The modes can then 
be spatially separated with a state-dependent potential 

m 

Finally, the state can be prepared with a 50/50 
beam splitter and a laser pulse corresponding to a state 
with low photon number variance. For obtaining (AiV) 2 
and (b''b'), a second beam splitter can be used, together 
with photon number measurements in the two modes. In 



oder to detect entanglement, assuming perfect destruc- 
tive interference at the second beam splitter, for the pho- 
ton source (AiV) 2 < N/4 — 7/8 is required. This can be 
obtained, for example, with a state with sub-Poissonian 
number statistics. 

Beside experimental considerations, the advantage of 
our approach is the ability to detect states in the vicinity 
of the entangled state J2J which is not detected by the 
method based on the correlation matrix [tj . The cor- 
relation matrix 7 contains the correlations of two pairs of 
conjugate single-party observables, which now we choose 
to be {R k } = {xa,Pa,x b ,Pb}- Here xa = (a + at)/y/2, 
PA — (a — at)/(v / 2i), and xb and pb are defined sim- 
ilarly for the b mode. For the state © the correlation 
matrix lkl = Tr{p(R k ~ (R k ))(Ri - (Jfc))} + Tr{p(R l - 
(Ri))(Rk — (Rk))} is obtained as 



' + 1 





N 








N + 1 
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N 





N+ 1 








N 





N + 



(19) 



The sufficient condition for inseparability is 7 — i J ^ 
where 7 is the correlation matrix corresponding to the 
partially transposed density matrix and J k i = i[Rk,Ri}- 
Here 7 — i J ^ thus the state is not detected as entan- 
gled. 

Moreover, with the simple method used for criterion 
||TJ described in the introduction, our criteria (@J, Q 
and l|14|) cannot be reduced to an entanglement witness. 
This is because they contain the variance of the particle 
number and (N) cannot be set to zero by single-party 
unitary operations. 



IV. CONCLUSIONS 

In summary, a simple inequality for the expectation 
values of observables was proposed for entanglement de- 
tection. Since only the measurement of easily acces- 
sible quantities (particle numbers and particle number 
variances) are needed, this approach may be feasible for 
detecting entanglement experimentally in Bose-Einstein 
condensates or in a two-mode photon field. 

Our method can be generalized for detecting other 
highly entangled states. First two operators must be 
identified which have the state as an eigenstate. Then 
a necessary condition for separability must be con- 
structed with the variances of these operators. Such a 
highly entangled state is for example the \N, 0) + |0, N) 
Schrodinger cat state which is the eigenstate of N and 
(a^b) N + (atf) N . 
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APPENDIX A: SINGLE MODE UNCERTAINTY 
RELATION 

1. Analytic calculation 

In this subsection we will prove Ineq. Ijllll . We will find 
a lower bound for the sum of the two variances (A p Na) 2 
and (A p a) 2 for any single-mode quantum state. This 
uncertainty relation is needed to find a lower bound for 
the sum of operator variances for two-mode separable 
states in Ineq. (jTJ). 

The first term on the left hand side of Ineq. (|llfl is 
zero for number states. The second term is zero for co- 
herent states. A/4 and a have a common eigenvector: for 
the state |0) the variances of both are zero. In order to 
find a non-trivial relation, the lower bound for the sum 
of the two variances must have at least one parameter. 
We choose this parameter to be (Na)„- For (Na) p > 
the operators Na and a do not have common eigenvec- 
tors and the sum of the two variances are bounded from 
below. 

The proof of Ineq. (|llf> is based on finding two lower 
bounds for the left hand side of Ineq. and then 

combining them. Let us denote (a) p = ^Ja (Na) „e 

where < a < 1. The first lower bound is 



where 



itfi 



R p > (A p a) 2 = (1 - a) (N A ) p =: B^Na, a). (Al) 
The second bound [Ig is obtained from 



(A P N A ) 2 + \[{A p x a ) 2 + {A pPa ) 2 ]-\ 

1 



R/j 



> ^2(A p N A ) 2 [(A p x A ) 2 + (A pPA ) 2 } - -. (A2) 

Here for the inequality X 2 + Y 2 > 2XY was ap- 
plied. Now, using the facts that (A p Na) 2 (A p xa) 2 > 
I (p A ) p 174 and (A p N A ) 2 (A pPA ) 2 > \ (x A ) p | 2 /4 we ob- 
tain 




y/a(N A ) p -±=:B 2 (N A ,a). 



(A3) 



From Ineqs. (|A1|) and i|A3() one can derive a higher lower 
bound for Ineq. by taking the maximum of these 

two bounds. It can be shown that 



B(N, a) := max [B 1 ,B 2 ] 



B±(N,a) if a > ctL, 
B 2 (N,a) if a< a L , 

(A4) 



1 



2VN 



(A5) 



Here oll is always non-negative, however, for small par- 
ticle numbers N < 1/4 it is larger than 1. 

The lower bound for Ineq. Ijlll) will be constructed 
by minimizing B(N,a) with respect to a. After some 
algebra one obtains 



min B(N, a) 



N- 



1 







if AT > 
if N < 



(A6) 



As stated in Sec. II. A, in order to use this result in the 
two-mode separability problem the bound should fulfill 
two criteria (i) L(N)+N 2 should be concave; (ii) L(Ni) + 
L(N 2 ) > L{Ni + N 2 ) + L(0). Equ. (HI does not fulfill 
(ii), thus a weaker bound satisfying this condition has 
to be chosen. L(N), as defined in Equ. p2|l . is such a 
bound. It coincides with the bound (|A6I) for A^ > 1/4 
while for A^ < 1/4 it is negative. 



2. Numerical verification 

In this subsection we prove by numerical calculations 
that Equ. (|12l) is a tight lower bound for Ineq. 1)110. We 
will determine the state vector minimizing the left hand 
side of Ineq. i|ll|) and the corresponding minimum with 
the constraint (a^a) = Na- 

The wave function is given in the number basis as 

|*) = £>|fc>. (A7) 

k 

The left hand side of Ineq. Ill 1(1 can be rewritten as 

R P = {EW^ 2 -^} 

+ [ n a- lE^+i^+if}- ( A8 ) 



Lagrange multipliers must be added in order to constrain 
the particle number to Na and keep the norm 1 

0({cm}> {C}>^i>M2) = R P - Vi{N A -^2\c k \ 2 k) 



/'2 



Cfc 



(A9) 



When minimized, all the derivatives of the function (|A9|) 
must be zero. Moreover, since R p ({ci~}) > R p ({\cf.\}) 
we can restrict our search for the minimum for real c^'s. 
Hence one obtains 



Cn+l = \ — = )Cn - [ \j— — )C„-1, (A10) 



71 + 1 
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FIG. 2: (diamonds) the minimum of R p (for a definition see 
Equ. Illl l if the particle number is constrained to Na- (solid) 
the analytic lower bound L(Na) defined in Equ. (1 1 2t . The 
inset shows the wave function in the number basis (diamonds) 
corresponding to the minimum of uncertainties for Na = 20. 
The results fit very well a Gaussian curve (solid). 



where A = (a) = J2 c kCk+iVk 
c„_i is not present for n 



j- 1 and the term with 
0. Based on Equ. IjAlOJI . 
from A, fi\ and \i 2 the unnormalized wave function can 
be constructed by setting cq = 1. As can be seen in 
Fig. |21 L{N A ) defined in Equ. ifPT)) is very close to the 
minimum found numerically thus it is a tight bound. The 
wave function minimizing R p is shown in the inset. In 
the number basis it fits very well a Gaussian curve even 
for small particle numbers. 



3. Generalized single mode uncertainty relation 

In this subsection we will prove Ineq. i|16fl . For w = 
the state minimizing R PiW is a coherent state, for w = 1 it 



is a number state. For intermediate w's the wave function 
giving the minimum interpolates between these two. The 
L w ((Na) p ) bound can be obtained, after inserting w and 
(1 — w) in the expression to be minimized, by following 
the same steps as in Appendix A.l. The two bounds 
found will be 



B ltW (N A ,a) = {l-w)(l-a)N A , 



B 2 , w (N A ,a) = \/w(l - w)aN A - 



1 — w 



(AH) 



The maximum of these two, B w (N,a), can be obtained 
knowing that Bi jW (N, a) > B 2 , W {N, a) if a > oil where 



v V 4JV(1 - w 

Hence the lower bound for R 



1 



AN(l-w) 
is obtained as 



(A12) 



min B,, 



yJw(l-w)(N 



i ifA>A L , 
if N < N L , 

(A13) 



where Nl = (1 — w)/Aw. 

As stated in Sec. II. B, in order to use these results 
in the two-mode separability problem the bound should 
fulfill two criteria (i) L(N) + wN 2 should be concave; 
(ii) L(Ni) + L(N 2 ) > L(Ni + N 2 ) + L(0). Equ. (lAT3ll 
does not fulfill (ii), thus a weaker bound satisfying both 
conditions has to be chosen. L W (N), as defined in l|17|l is 
such a bound. It coincides with Equ. (|A13|I for N > Nl 
while for < Nl it is a linear function of N and it 
is negative. The function giving L W (N) for N > Nl 
(top line in Equ. (|A13p ) cannot simply be extended to 
N < Nl as it was done in Appendix A.l for the simpler 
uncertainty relation, since in this case (i) would be not 
satisfied. 
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